Fermi surface pockets in ortho-II YBa2Cu306.5: the origin of quantum oscillations? 
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In this paper we explore whether the quantum osciUation signals recently observed in ortho- 
II YBa2Cu306.5 (OII-Y123) may be explained by conventional density functional band-structure 
theory. Our calculations show that the Fermi surface of OII-Y123 is extremely sensitive to small 
shifts in the relative positions of the bands. With rigid band shifts of around ±30 meV small tubular 
pockets of Fermi surface develop around the Y point in the Brillouin zone. The cross-sectional areas 
and band masses of the quantum oscillatory orbits on these pockets are close to those observed. 
The differences between the band-structure of OII-Y123 and YBa2Cu408 (Y124) are discussed with 
reference to the very recent observation of quantum oscillations in Y124. 



The nature of the normal state of the high temper- 
ature cuprate superconductors (HTC) has been a topic 
of intense discussion ever since their discovery. The un- 
usual temperature dependence shown by the resistivity 
and Hall coefficient for example, and how these evolve 
as a function of doping, has led to a wide range of ex- 
otic theories of the nature of the normal state.-l Many 
of these theories depart substantially from the standard 
Fermi liquid model of a metal, particularly in the under- 
doped region of the phase diagram. 

The recent report^ of the observation of Shubnikov- 
de Haas (SdH) oscillations in the Hall and longitudinal 
resistivities of underdoped ortho-H YBa2Cu306.5 (OII- 
Y123) was a surprising and potentially extremely impor- 
tant result. These oscillations suggest that even in the 
underdoped region some well-defined pockets of Fermi 
surface exist. The observed orbits have a frequency of 
~ 530 ± 20 T and a mass, m* — 1.9 ± O.lwe {rrie is the 
free electron mass). A single pocket with this dHvA fre- 
quency corresponds to a small (~ 2%) fraction of the Bril- 
louin zone. Identifying the origin of this orbit is clearly 
very important for its interpretation. In this paper, we 
discuss whether the SdH signals could come from small 
pockets of Fermi surface predicted by conventional den- 
sity functional band-structure calculations. 

Our understanding of the electronic structure of the 
cuprates has, to date, mostly been led by angle resolved 
photoemission spectroscopy (ARPES). The resolution of 
this technique has advanced rapidly over the past decade 
and had provided many key insights into the physics of 
HTC.^ More recently, the strong angle dependent mag- 
netoresistance of strongly overdoped Tl2Ba2Cu06+5 has 
been used to extract information about the Fermi surface 
and scattering rate.^^ The shape and size of the Fermi 
surfaces measured by both these techniques are gener- 
ally in very good agreement with conventional density 
functional theory (DFT) band-structure calculationsi^i^ 

Quantum oscillatory effects, e.g. the de Haas-van 
Alphen (dHvA) and Shubnikov-de Haas effects, are very 
powerful probes of the Fermi surface properties of a 
metal. Unlike ARPES, they probe the bulk of the ma- 
terial (and so are not sensitive to surface defects) and 
are a true three dimensional probe of the quasiparticles 
at the Fermi level. The analysis of the amplitude and 



frequency of the signals and how these change as a func- 
tion of magnetic field and temperature, gives information 
about the Fermi surface and the quasiparticle masses and 
scattering rates. The main constraint is that impurities 
severely attenuate the SdH/dHvA signal, according to 
the Dingle factor Rd ~ cxp(-^^g^) (here 7rfc|, ~ A the 
cross-sectional area of the dHvA orbit and ^ is the or- 
bitally averaged mean free path). The small size of the 
orbit in OII-Y123 was a key factor in it being observable. 

Band-structure calculations of fully oxygenated ortho- 
I YBa2Cu307 (OI-Y123) show a small tubular pocket 
of Fermi surface near the S point which derives from 
the 'CuO chain bands'^"^!^ (see below). This could be 
a simple explanation for the origin of the observed SdH 
orbit. However the band-structure calculations of Bas- 
cones et alr^ indicate that this small pocket is absent 
in OII-Y123. Another intriguing possibility is that the 
SdH signal comes from pockets close to the nodal regions, 
which may have formed because of Fermi surface recon- 
struction. In this paper we calculate the band-structure 
of OII-Y123 using conventional density functional theory 
techniques and investigate the possible extremal orbits 
which could give rise to the observed orbit. We contrast 
this with similar calculations on the double-chain cuprate 
YBa2Cu408 (Y124) in which SdH oscillations have also 
recently been observed i^^i^^ 

Our calculations were carried out using the Wien2K 
packagCfi^ which is an implementation of a full-potential, 
augmented plane wave plus local orbital scheme. We 
used a generalized gradient approximation form for the 
exchange correlation potential, ^"''^^ and the crystal struc- 
ture of YBa2Cu306.5 determined by Grybos et al.i^^ 
A dense fc-mesh of ^ 10^ points in the full Brillouin zone 
(19 X 39 X 12, or 1400 points in the irreducible wedge of the 
Brillouin zone) was used for the self-consistency cycle. 

The calculated band-structure in the basal plane at 
the center of the Brillouin zone {k^ — 0) is shown in Fig. 
[1] The main features can be easily related to that of 
fully oxygenated OI-Y123 by band folding. In the first 
panel of the Fig. [T{a) we show our calculation of the 
band-structure of OI-Y123 which is almost identical to 
previous calculations (e.g., Refs.7,9). The 01 results have 
then been folded down into the smaller OH Brillouin zone 
produced by doubling the a lattice parameter (panel b). 
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FIG. 1: (Color online) Panels (a-c): The band-structure of YBa2Cu307_i. (a) OI-Y123, (b) OI-Y123 folded in to the reduced 
zone of OII-Y123 (c) OII-Y123. Panels (d-f) Fermi surfaces of YBaaCusOr-i . (d) OI-Y123, (e) OI-Y123 folded, (f) OII-Y123 
(with AEf^+20 meV). 



It can been seen that this procedure reproduces many 
of the features of the full OII-Y123 calculation (panel 
c). The effect of the band folding on the Fermi surface 
is also shown in the figure and again compared to the 
results from the OITY123 calculation (panels d-f). 

Moving along the X-S line the first band to cross 
the Fermi level is mainly due to the CuO chain (Fig. 
[T];). It can be seen that this is very close to being one- 
dimensional and has much less dispersion compared to 
the 01 result. Note that there is only one chain band 
crossing in the Oil calculation compared to two in the 
folded 01 calculation because Oil only has one conduct- 
ing CuO chain per unit cell. The next four crossings are 
due to the Cu02 planes. As in the 01 calculation there 
is a sizeable splitting between the bonding and antibond- 
ing Cu02 bands 210 meV along the X-S line at Ep). 
Each of these two bands are further split into two by 
the additional 2a periodicity. This is most evident close 
to the S point. We find that this sphtting is very small 
at Ep along the X-S line, being 17 ± 2 meV. In the 
OTY123 calculation a fairly flat CuO/BaO band, arising 
mainly from hopping between chain oxygen sites via the 
apical O site in the BaO layer below it, passes through 
Ep close to the S point and gives rise to a small tubular 
(quasi- two-dimensional) hole pocket Band folding 

moves this same pocket to the Y point in the small (Oil) 
zone. In the full OII-Y123 calculation, this band also 
crosses Ep and again forms a tubular hole pocket close 
to the Y point. In addition, a second band with the same 
character is very close in energy, ~ 35 meV lower at Y. 

Our results are similar to a calculation of OITY123 re- 
ported by Bascones et ali^ The main differences are that 
the splittings of the Cu02 bands due to the 2a period- 
icity are larger and near to Y the CuO/BaO bands are 
more closely spaced (and ^ 50 meV lower in energy) than 
in our calculation. Bascones et al. used a TB-LMTO- 
ASA method which may be less accurate for a non-close 



packed structure like OII-Y123 than the full potential 
method used here. In any case, these differences would 
not change our essential conclusions. 

The shape and size of the two plane sheets predicted 
for OTY123 have been verified by ARPES measurements 
on YBa2Cu3 06.95.^^ There is also some weak evidence for 
the quasi- ID chain sheet but there is no sign of the hole 
pocket centered on S. However, it should be mentioned 
that ARPES measurements on YBCO are complicated 
by a surface feature originating from the CuO chains, 
which may limit the resolution for Fermi surface sheets 
with significant chain character 

Although DFT calculations usually correctly describe 
the general features of the band-structure, there are of- 
ten small discrepancies with respect to the relative posi- 
tions of the bands when compared to experiment. This 
is true even if the calculation has been converged to the 
meV level. Sr2Ru04 is an example of an oxide mate- 
rial, with a similar structure to the cuprate supercon- 
ductor La2_2:Sr2:Cu04, which has been studied in detail 
by dHvA measurements. To get agreement with band- 
structure calculations the bands need to be shifted by 
^^40 meV, in opposite directionsi^ Even in the rela- 
tively simple material MgB2, shifts of order 100 meV are 
needed 

In Fig. [2] we show how the Fermi surface of OIT 
Y123 changes as the Ep is varied by ±35 meV. These 
small rigid band shifts correspond to adding (remov- 
ing) 0.05 (0.07) electrons per Cu atom, although a sim- 
ilar result could be obtained without doping by moving 
the plane and chain bands in opposite directions. For 
A£^i?=+35mcV the Fermi surface consists of two large 
hole-like tubular Cu02 sheets centered on S, plus three 
quasi-one-dimensional sheets (one from the chains and 
two from the planes). As Ep is reduced a small hole-like 
pocket develops near the Y point, the origin of which (as 
discussed above) is the same as that of the pocket near 
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FIG. 2: (Color online) Evolution of the Fermi surface of OII-YBa2Cu306.5 with Fermi level shift AEf- The figure shows 
two-dimensional cuts within the basal plane {k^ =0), as indicated by the symmetry labels. The primary quantum oscillation 
orbits {Fn) are marked on the +20 meV and -20 meV panels. 



the S point in OI-Y123. Further reduction of Ep re- 
sults in this pocket growing in size and then merges with 
the antibonding Cu02 plane sheet. As Ep is further 
reduced the second CuO/BaO band passes through the 
Fermi level giving rise to another pocket. Eventually this 
merges with the bonding Cu02 plane sheet. Both of the 
pockets that open close to Y are quasi-two-dimensional 
(tubular) sheets, with relatively weak warping along the 
c direction (see Fig. [3]). Similar results to these have 
recently been report by Elfimov et al.^ 

The above illustrates the extreme sensitivity of the 
Fermi surface to the relative positions of the bands and 
means that subtle changes to the doping could result in 
small Fermi surface pockets being formed. It is possible 
that one such pocket could give rise to the quantum os- 
cillation signals observed by Doiron-Leyraud et al.^ To 
investigate this possibility more fully we have calculated 
the quantum oscillation frequencies [F = (11/2^6) A] and 
band masses [mi, — ^%^]- The frequencies shown in the 
bottom panel of Fig. 2| (F3 and F4) are from the small 
hole pockets discussed above, whereas those in the up- 
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FIG. 3: (Color online) Three dimensional image of the cal- 
culated Fermi surface of OITY123 showing the small pockets 
centered on the Y points for the two different Fermi level 
shifts indicated. The shading indicates the band character 
C=Chain like (including apical oxygen), P=plane like. 
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FIG. 4: (Color online) Extremal dHvA frequencies versus 
Fermi level shift (solid symbols). The corresponding band 
masses are shown with open symbols (right-hand axis). 



per panel (F\ and F-i) are from the main Cu02 sheet 
surfaces. Each of the tubular sections of Fermi surface 
has a minimum and maximum extremal area which have 
been labeled i^" and respectively in the figure. As 
expected, the extremal areas (dHvA frequencies) of both 
hole pockets vary strongly with the Fermi level shift and 
have maximum frequencies of 400-600 T. The band mass 
of the pocket F3 is roughly constant as a function of 
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FIG. 5: (Color online) Band-structure of YBa2Cu408 (Y124). 
The symmetry labels appropriate for the Ammm space group 
of Y124 are the same as in Ref|2^. The path FZJI is approxi- 
mately equivalent to the FXSY path used in the other figures 
for the Pmmm space group of OI/II-Y123. 

AEp with an average value of ~ 1.5 mg. The mass of 
the other pocket orbit, F4 is much more variable, rang- 
ing from ~ 1.5 for large AEp to ^ 3 as AEp approaches 
zero. A similar calculation for OI-Y123 was reported in 
Reftii. 

The experimentally observed orbit frequency {F = 
530 ± 20 T) is at the upper limit of the calculated range 
for both F3 and F4. The observed mass of the orbit 
(to* = 1.9 ± O.lTOe) seems to be more compatible with 
F3 than F4. We expect that the observed quasiparticle 
mass should be enhanced compared to the band mass 
[m* = (1 -|- X)mb\ because of many body interactions 
(electron-phonon and electron-electron) which are not in- 
cluded in the present calculations. Although it is ex- 
pected that A will be relatively large (~ 1 — 2) on the 
Cu02 planes^^ it could be substantially less on these 



sheets of Fermi surface without significant plane charac- 
ter. Hence, the calculated mass of F3 could be consistent 
with experiment. Experimental studies of the pressure or 
doping dependence of the SdH frequencies in Y123 may 
help determine this interpretation is correct. An alterna- 
tive approach is to study other cuprate materials where 
the CuO/BaO band is either absent or is at markedly 
different energies. 

Very recently two groups^»i^ have reported the ob- 
servation of SdH oscillation in Y124 with a compara- 
ble frequency to those in OII-Y123. Y124 is an intrin- 
sically underdoped cuprate with a similar structure to 
OI-Y123, but with a double chain layer and a doubled 
unit cell along the c-direction. Our calculations of the 
band-structure, which are essentially the same as re- 
ported previously r^i^ show that the CuO/BaO band 
which gives rise to the pockets in Y123 is now around 
400 meV below Ep at the zone corner J (marked with 
an asterisk in Fig. [5]). Hence, it is very unlikely that this 
band could be responsible for the experimental observa- 
tions. In fact, we not find any orbits within AEp = ±400 
meV which have frequencies close to experiment. 

In summary, we have presented calculations of the elec- 
tronic structure of OII-Y123, and have discussed possi- 
ble origins of low frequency quantum oscillation signals. 
Two possible orbits, both arising from chain/apical oxy- 
gen bands have been identified, which with a small shift 
in the relative energies of the bands, have extremal ar- 
eas compatible with the observed SdH signals. The band 
masses of the quasiparticles on at least one of these or- 
bits is roughly consistent with those observed. It seems 
unlikely however that the reported orbits in Y124 can be 
explained by a similar band-structure approach. 

We thank S.B. Dugdale, N.E. Hussey, and S.M. Hay- 
den for useful comments. 
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